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ABSTRACT 

The orbit of a planet is surrounded by a chaotic zone wherein nearby particles' orbits are 
chaotic and unstable. Wisdom (1980) showed that the chaos is driven by the overlap of mean 
motion resonances which occurs within a distance ((5a/a)chaos ~ of the planet's orbit. 

However, the width of mean motion resonances grows with the particles' eccentricity, which 
will increase the width of the chaotic zone at higher eccentricities. Here we investigate the 
width of the chaotic zone using the iterated encounter map and N-body integrations. We find 
that the classical prescription of Wisdom works well for particles on low-eccentricity orbits. 
However, above a critical eccentricity, dependent upon the mass of the planet, the width of the 
chaotic zone increases with eccentricity. An extension of Wisdom's analytical arguments then 
shows that, above the critical eccentricity, the chaotic zone width is given by {Sa/a)chaos ~ 
1.8ei/V^^^ which agrees well with the encounter map results. The critical eccentricity is 
given by Ccrit ~ This extended chaotic zone results in a larger cleared region 

when a planet sculpts the inner edge of a debris disc comprised of eccentric planetesimals. 
Hence, the planet mass estimated from the classical chaotic zone may be erroneous. We apply 
this result to the HR 8799 system, showing that the masses of HR 8799 b inferred from the 
truncation of the disc may vary by up to 50% depending on the disc particles' eccentricities. 
With a disc edge at 90 AU, the necessary mass of planet b to cause the truncation would be 
8-10 Jovian masses if the disc particles have low eccentricities (< 0.02), but only 4-8 Jovian 
masses if the disc particles have higher eccentricities. Our result also has implications for the 
ability of a planet to feed material into an inner system, a process which may explain metal 
pollution in White Dwarf atmospheres. 

Key words: chaos — celestial mechanics — planets and satellites: dynamical evolution and 
stability — circumstellar matter — stars: individual: HR 8799 



1 INTRODUCTION 

Planets strongly affect the dynamics of bodies on nearby orbits. 
Bodies coming within around one Hill's radius are strongly per- 
turbed on the synodic time-scale and can either be scattered, result- 
ing eventually in collision with the star or ejection from the sys- 
tem, or collide with the planet. However, beyond this region there 
is a chaotic zone of unstable orbits. Here, orbits diffuse through 
phase space and may eventually find themselves intersecting the 
planet's orbit, which will again result in a strong scattering or colli- 
sion event. This unstable chaotic zone is thought to be responsible 
for sculptin g the inner edges of debris discs such as that aroun d 
Fomalhaut jouilleij|2006ah iKalas et al]|2008l : IChiang et alj|2009h . 
with the location and shape of the disc edge being determined by 
the mass and location of the planet. The zone may also provide a 
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reservoir for feeding bodies from an outer disc into the inner sys- 
tem. This is likely to be important in the late stages of stellar evolu- 
tion, when scattering of planetesimals by planets has been invoked 
to explain metal pollution observed in some White Dwarf atmo- 
spheres as well as hot discs around them (Z uckerman et al. 20031; 
iFarihi. Jura & Zuckermarj|2009l ; lBonsor. Mustill & Wvatll201 ll) . 

Chaos in dynamical s y stems is often d riven by overlapping 
resonances i lChirikov|[T97^ . IWisdod Jl980h derived an analytical 
formula for the extent of a planet's chaotic zone based on overlap- 
ping first-order mean motion resonances, finding 



[Sa/a] 



chaos C/i 



2/7 



(1) 



where /i = rripi/m* is the ratio of planetary to stellar mass and 
the coefficient c = 1.3. This has been co nfirmed numerically , 
both using iterations of the encounter map touncanetai]|l989l 
see §[2]below) and in full N-body integrations (Chian g et alj200^ . 
although the coefficient c when determined numerically is some- 
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what higher. However, the result of IWisdon] ( ll98Cll) is only valid 
for low-eccentricity particles, since the widths of mean motion res - 
onan ces increase at larger eccentricities j Murrav & DermottI 19991) . 
lOuille n & Faber ( 2006) extended Wisdom's result to consider the 
behaviour of particles orbiting close to an eccentric planet with low 
free eccentricities, which is relevant for modelling the Fomalhaut 
disc. They showed that the width of the zone for particles orbit- 
ing close to an eccentric planet with low free eccentricities is the 
same as that for low-eccentricity particles orbiting close to a circu- 
lar planet. 

The effect the particles' free eccentricities have on the width 
of the chaotic zone has not been well explored. Since the Kuiper 
Belt Objects of our own Solar System can have free eccentricities 
of 0. 1 and above, one might expect the debris discs of some extra- 
Solar systems to be stirred to a similar degree. Isonsor et al] ^201 ih 
found, in N-body integrations, that the width of the chaotic zone 
for more eccentric particles (e up to 0.1) appeared to be greater 
than that for less eccentric particles. The small number of particles 
in the N-body integrations did not permit definite conclusions to 
be drawn, but they confirmed the increasing width using the com- 
putationally more efficien t iterated encounter map ( Du ncan et al.l 
1 19891) . I Bonsor et al. l j20H speculated that the increasing width of 
the zone were due to the particles' pericentres entering the chaotic 
zone. This however seems unlikely since the chaos is caused by 
resonance overlap and that depends on the particles' semi-major 
axes and not directly on their pericentres. Here we demonstrate that 
the increasing chaotic zone width is instead due to the increasing 
width of mean motion r esona nces as particles' ecc entricities are in- 
creased. lOuilleij ( I2OO7I) and lChiangetal] j2009l) investigated the 
effects of non-gravitational effects such as collisions and radiation 
forces on the inner edge of the disc, and here we investigate the 
effects of the particles' eccentricity. 

One reason why it is important to characterise this chaotic 
zone properly is that currently much interest revolves around the 
effects that perturbing planets can have on debris disc morphology, 
particularly with a view to determining or refining the parameters 
of known or suspected planets. Accordingly, the /i^'''^ law has been 
used to constrain the mass and semi-major axis of planets that may 
be sc ulpting the inner edges of debris discs {e.g., iMoerchen et al.l 
I2OI i ll), to verify the consistency o f planet and dis c parameters deter- 
mined through other means (e.g., ISu etalj2009|K and to aid setting 
up and interpreting N-body integrations (e.g., Chiang et al. 2009) 
studying the sculpting of debris disc edges. Hence, refinements to 
and deviations from this law are of prime importance for under- 
standing the interaction of planets and debris discs. Given the large 

uncertainties which attend determinations of masses of planets de- 
i 1' — 

tected by direct imaging cKalas et al . 2008), dynamical constraints 
are very valuable and use of accurate dynamical models should be 
made. In §|4j2]of this paper, we discuss the HR 8799 system and 
how mass estimates of planet b, thought to be sculpting the inner 
edge of the cold debris disc, are affected by deviations from the 
^i^/^ law. 

This paper is organised as follows. In §|2]we use the encounter 
map and N-body integrations to explore the chaotic zone as a func- 
tion of particle eccentricity. In §|3]we derive the width of the chaotic 
zone for eccentric particles using a resonance overlap criterion. In 
§|4]we discuss the implications of this for studies of debris discs 
sculpted by planets, particularly the HR 8799 system, and plane- 
tary systems orbiting White Dwarfs. 



2 NUMERICAL INVESTIGATION 

We investigated the structu re of the chaot i c zone using the iterated 
encounter map described in lDuncan et T his uses the ap- 

proximate solution to Hill's Equations calculated bv lHenon & PetitI 
( 1986), which assumes a moderate separation between the planet 
and particle, to calculate the impulsive change to a particle's orbital 
elements due to a conjunction with the planet every synodic period. 
We determined whether orbits were regular or chaotic by examin- 
ing the Fourier Transform of the eccentricity evolution. Regular 
orbits have power spectra with a few well-defined peaks, while the 
more chaotic an orbit the more peaks its power spectrum displays 
(this is a consequence of the phase space of a chaotic dynamical 
system being densely filled with periodic orbits). We chose a crit- 
ical number of peaks, a peak being defined as where the slope of 
the power spectrum changes from positive to negative, to decide be- 
tween regular and chaotic trajectories. For the investigations below, 
with 10* iterations, we found that 1000 peaks distinguished regular 
and chaotic orbits over the whole range of parameters considered. 
This class ification method is similar to the spectral number chaos 
indicator jZhou et al .11201 ll) . 

In this way we then investigated the nature of trajectories as 
a function of the particles' semi-major axis e = (a — api)/api 
and eccentricity e, and the planet:star mass ratio /i. The results are 
shown in Figure[T] In these plots, each point in the grid represents 
100 particles, with initially random values of longitude of pericen- 
tre and mean anomaly. The greyscale shows the fraction of particles 
whose orbits were classed as chaotic, with black representing 100% 
regular and white 100% chaotic. In the region at the top left comer 
above the black line the particles are on planet-crossing orbits and 
hence the encounter map is invalid. These particles were classed as 
chaotic since they will be unstable anyway except perhaps for a few 
protected from close encounters by resonances. 

From Figure [T] we see that the edge of the chaotic zone is 
a function of eccentricity, with the more eccentric particles being 
chaotic at greater separations from the planet than their less ec- 
centric brethren. There are two regimes: a low-eccentricity regime 
in which the chaotic zone is independent of eccentricity, and a 
moderate-eccentricity regime where the width of the chaotic zone 
increases with eccentricity. The former regime is that identified by 
Iwisdod ( fl980h . The critical eccentricity separating these regimes 
increases with /x. The width at low eccentricity is in go od agree- 
ment with the result from the resonance overlap criterion l lwisdomi 
Il980l shown by a vertical blue line), although this systematically 
under estimates it, as previous studies have shown (Duncan etalj 
Il989l) . The increased chaotic zone width at larger eccentricities is 
likely due to the increasing width of the mean motion resonances 
at higher eccentricities, which can clearly be seen as intrusions of 
chaotic motion into the regions of regular motion. This increasing 
width can be derived analytically, as shown below. 

To verify the accuracy of the encounter map results, we ran an 
N-body integration using the hybrid algorithm included in the Mer- 
cury package (Chambers 1999). 961 test particles were placed on 
orbits exterior to a planet of mass ^ = 10^* and their orbits inte- 
grated for 10 Myr. These parameters correspond to Panel (f) of Fig- 
ure [T] Their final eccentricities and semi-major axes are shown in 
Figure|2] Initially they fill the shaded region, but many are removed 
from regions closer to the planet and scattered onto orbits which 
almost intersect the planet's (shown by the dashed black line). The 
vertical solid b lue line shows the chaotic zone width according to 
IWisdodfl980l) . More eccentric particles are destabilised at greater 
semi-major axes than less eccentric ones, as was seen in the results 
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Figure 1. (a)-(f): Extent of chaotic zone as a function of particle eccentricity and semi-major axis and planet mass. The right-hand side of each plot is at two 
Hill radii from the planet. The vertical blue line shows the extent of the classical chaotic zone (Wisdom, 1980). Above the dashed black line, the particles 
ai'e on planet-crossing orbits; these pailicles were all classed as chaotic regardless of the map behaviour since they will experience close encounters with the 
planet. The solid red line shows the analytical estimate of the chaotic zone width as a function of eccentricity (Equation llO) . 



from the encounter map. Hence, the encounter map appears to be 
giving an accurate picture of the dynamics. Additionally, the region 
above the dashed black line is almost totally devoid of particles, due 
to scattering by the planet (these particles have either collided with 
the planet or are in the unstable population just below the black 
line which is in the process of being scattered). This justifies our 
labelling of planet-crossing orbits as "chaotic" in Figure[T] 



3 ANALYTICAL DERIVATION 

The standard derivation of the chaotic zone width JWisdonj[T980l : 
lOuillen & Fabeill2006 ) takes the widths of the mean motion reso- 
nances to be independent of eccentricity. Working from the reso- 
nance Hamiltonian 

+ I3J - j'^^ cose, (2) 

where J cx is the canonical momentum, 6 is the resonant argu- 
ment, and /3 is a parameter measurin g distance to the nom inal reso- 
nance location (see lOuille n 2006b; Mustill & Wvattl201 ih . one can 
show that the width of a first-order resonance at low eccentricity is 
approximately 

iSa/a)rcs = 3.2^1^^'' f^'^ (3) 

dWisdomll 19801) . Here j numbers the resonance, j : (j — 1). In the 
limit of large j, adjacent resonances are separated by a distance 

(Sa/aUp = 2/{3f). (4) 
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Figure 2. Fates of 96 1 particles after after a 1 Myr N-body integration. Ini- 
tially they are distributed over the grey rectangle and their final semi-major 
axes and eccentricities are plotted as crosses. Also shown are the same lines 
as in Figure 1: Dashed black; planet-crossing orbits; Solid blue; classical 
chaotic zone; Solid red; extended chaotic zone. Dynamical clearing above 
the red line is apparent. 
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Figure 3. Black, solid line and crosses: The chaotic zone width at a mod- 
erate eccentricity of e = 0.01 as a function of planet:star mass ratio. Points 
show the numerical results from the grids in Figure [T] and the line shows 
the regression line, given by (<5ii/ii)chaos = 0.86/i"'^'^*, close to the an- 
alytical result predicted by Equation[To] Blue, dotted line and stars: The 
chaotic zone width at zero eccentricity. Here the regression line (dotted) 
has a slope of 0.278 ± 0.0077, close to the 2/7 predicted by the analytical 
arguments of Wisdom (1980). 



Equating these two distances gives tlie outermost resonance where 
overlap occurs as 



j = 0.51m 



-2/7 



(5) 



Since a — Opi = 2/{3j) for large j, the chaotic zone width is thus 
given by 



(5a/a) 



l.S^i 



2/7 



(6) 



This scaling is succ essfully reproduced by numerical iterations of 
the encounter map ( Duncan et aljl98^ . and Figure[3]below). 

However, the resonance width grows with particle eccentric- 
ity, as can be seen in panel (a) of Figure [T] This means that res- 
onances that do not overlap at low eccentricities may overlap at 
higher eccentricities, and cause the chaotic zone to be wider. When 
eccentricities are not small, the width of the resonance can be 
found by considering the width of the resonant libration region af- 
ter the resonant separatrix forms, which is given by SJ = 4(4 J) 
jMurrav & Dermotlll 19991) in terms of the canonical momentum J. 
Now J is related to e by 



J = kj{mpi/m(s) ^^^{m^,/m@f^^e^, (7) 

where fc, ^ tj./mp 

dMustiU & Wvattl201 iklouillen & FabeJlOOd) . and so 

5J = l.S^-^/^i'/^eJe 



7 (3,7'mp/2mffi)'^^^/2 in the limit of large j 

(8) 



Using the relationship between eccentricity cha nge and semi-major 
axis change in a resonance Sa — 2jae5e JMurrav & Dermotd 
Il999t) . we then have 



for the maximum libration width in semi-major axis. Equating this 
to the resonance separation and redoing the previous derivation 
then gives 



(10) 



as the width of the chaotic zone at higher eccentricities. Note the 
two implications of this equation: first, the chaotic zone width 
grows weakly with eccentricity. The dependence is shown in Fig- 
ure [T] as solid red lines. The agreement with the results from the 
encounter map is striking. We also show the dependence on Fig- 
ure (2] where we can see many unstable particles above the solid 
red line. Second, the fi^^'^ dependence of the chaotic zone width 
on planet mass, valid in the low-eccentricity regime, is replaced 
with a fj,^^^ dependence. To verify this, in Figure [3] we plot the 
width of the chaotic zone at an eccentricity of e — 0.01 as a func- 
tion of fi. The width was determined by finding the innermost cell 
where all trajectories were regular from the grids presented in Fig- 
ure [T] The regression line through these points gives a dependence 
((Sa/a)chaos = 0.86/^"'^'^*, in excellent agreement with the ana- 
lytical result from Equation 11 01 Thus, for higher eccentricities, the 
chaotic zone width is given by a jj,^''^ scaling, not the more famil- 
iar /i^''^. For comparison, the width of the chaotic zone at e = 
as a function of /i is shown on the same Figure. The regression line 
here gives ((5a/a)chaos = 1.35f/''^^^ , in excellent agreement with 
Wisdom's original result. 

Equation[TO]also allows the critical eccentricity separating the 
two regimes to be estimated. Equating ((5a/a) chaos obtained from 
Equation[TOl for eccentric particles, with the value from Equation|6] 
for low-eccentricity particles, gives 



Gcrit 



: 0.21^I 



3/7 



(11) 



growing with ^, as is seen in Figure [T] Note that e = 0.01 is the 
critical eccentricity for /i — 10^^, and indeed the widths of the 
classical and extended zones coincide at this point (Figure[3). 

Note that there is a maximum eccentricity for which this 
derivation of the extended chaotic zone is valid. This occurs where 
the solid red line in Figure[T]marking the boundary of the extended 
chaotic zone intersects the dashed black line denoting the region of 
planet-crossing orbits. At eccentricities above this, particles will be 
removed directly by close encounters rather than via chaotic diffu- 
sion. The maximum eccentricity is given by 

e - 2 lu^/"* 



and in this regime particles are unstable if 

Sa/a < e. 



(12) 



(13) 



(Sa/a) 



7.78e 
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(9) 



In the next Section we adopt the slightly more conservative con- 
dition that particles are also unstable if they approach within one 
Hill's radius of the planet. 



4 DISCUSSION 

4.1 The shape of debris disc inner edges 

This result has important implications for the interactions of plan- 
ets and debris discs. Discs showing inner clearings are commonly 
explained by the existence of a planet removing particles from 
its chaotic zone (e.g., IOuillenii20()6al) . The mass and semi-major 
axis of the pla net determine both the l ocation and the sharpness of 
the in ner edge dOuillen & Fabej2006l : IOuillerj2006iJ:IChiang et all 
I2OO9I) . Since our investigations above have revealed that the extent 
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Figure 4. Surface densities at tlie inner edges of debris discs cleared by a 
planet at 1 AU. Top: Comparison of the encounter map results (solid) with 
the prescription of Wisdom (1980; dotted). Vertical dotted lines show the 
location of the classical zone edge in semi-major axis; irregular dotted lines 
show the planetesimal distribution when their eccentricities are accounted 
for. Maximum particle eccentricities are 0.1. Planet masses are /i = 10~^ 
(black), pL = 10^* (blue) and = 10~^ (red). Bottom: Disc profiles 
according to the encounter map results for different planet masses and plan- 
etesimal eccentricities. Planet masses are as above; maximum particle ec- 
centricities are 0.01 (solid), 0.03 (dashed) and 0.1 (dotted). 



of the chaotic zone differs from the classical prescription, we now 
explore its ef fect on the nature o f inner edges to debris discs. As 
discussed bv iChiang et alj j2009h . to correctly model the shape of 
a disc edge it is necessary to account not only for particles located 
in semi-major axis at the edge of the chaotic zone, but also ones 
further out whose eccentricities may cause them to affect the radial 
distribution of particles at the disc edge. Accordingly, we generate 
surface density profiles for discs in the following way. We consider 
simple model discs where 5 x 10* planetesimals are uniformly dis- 
tributed in eccentricity between and a maximum em and semi- 



major axis between 1 AU and 2 AU. Particles are removed if their 
orbits under iterations of the encounter map are chaotic, or if their 
periapsis comes within one Hill's radius of the planet. We then sam- 
ple 100 mean anomalies from each particle's orbit, and sum the 
particles in radial bins to find the surface density profile. We also 
vary the particles' eccentricities and semi-major axes according to 
the encounter map results, because the particles do not remain at 
their initial a and e even if their orbits are regular. 

The results, for mass ratios /i — 10^'^, 10^* and 10^^ and 
a maximum eccentricity of em = 0.1, are shown in the top panel 
of Figure |4] as solid lines. We see that the shape of the inner disc 
edge depends on both the planet's mass, with more massive planets 
making the edge shallower. We show on the sam e plot with dotte d 
lines the results from the original prescription of Iwisdonil ( Il98(]h . 
showing a discrepancy both in the location of the edge and in the 
shape of it. The discrepancy is greater for larger planet mass. In 
the bottom panel of Figure |4l we show the results from the en- 
counter map iterations for a range of planet masses and particle ec- 
centricities. More massive planets and more eccentric particles lead 
to shallower slopes. In addition, fluctuations in the surface density 
outside the cleared zone caused by mean motion resonances can be 
seen (see also Ouillen 2007), and these are more apparent for lower 
eccentricities. The shallowness of the edge means that its precise 
location is ill defined, but if it is attempted to be quantified by a 
single statistic, for example the half-maximum, then it should be 
noted that this value depends on the particles' eccentricities, lead- 
ing to a degeneracy between planet mass and particle eccentricity. 

In this subsection, and the subsequent, we consider only the 
gravitational influence of the planet on the planetesimals which 
the disc comprises. It is likely that in many debris discs the ec- 
centricities of planetesimals are high, if they h ave been dynam- 
ically excited by planets' secular perturbations ^Mustill & Wvatll 
1 20091) or by large plane tary embryos embedded in the disc 
I Kenvon & BromIevl20ICI) . In such a disc collisions between plan- 
etesimals will be destructive due to the high relative velocities. In 
such a disc eccentricities will remain high as ongoing stirring pro- 
cesse s will dominate over co Uisional damping (Kenvon & BromlevI 
|2004| ; [Shannon & Wuir20I ih . We do not e that there may be discs 
where the relative velocities remain low faeng & Tremainell20IOl) 
in which eccentricities may be continually damped by collisions. 
The truncation of such discs will be governed by both dynamical 
clearing from the classical chaotic zon e and viscous spreading of 
the disc due to the inelastic collisions jOuilleij|2007h . We do not 
consider such discs here. 



4.2 Case study: HR 8799 

As an example of the effects of the extended chaotic zone, consider 
the HR 8799 system. The 1.5 Solar mass star is orbited by 4 known 
planets and at least two debris belts. ISu et al. 1 12009) estimated that 
the inner edge of the outer disc lies at 90 AU based on fitting the 
disc's spectral energy distribution (SED); the inner edge of the disc 
is not resolved. It is thought that the outermos t planet, HR 8 799 b, 
at a projected distance of 68 AU f rom the s t ar jMarois et alfcoOSh , 
is responsible for the truncation. ISu et al.l j2009l) argued that the 
radius of the truncation, the location of planet b at 68 AU, and a 
planet mass of 10 Jovian masses, at the upper end of the range sug- 
gested by photometry (Marois et al. 2008), are all consistent with a 
scenario in wh ich the tr uncation is due to planet b's chaotic zone. 
Ouantitativelv. lSu et alj(l2009i) took t he chaotic zone to be given by 
1.4/i'^^^ as given by iMalhotral ( [T998h . 

As we have seen, this relationship underestimates the width 
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of the chaotic zone at high eccentricities. If the disc particles' ec- 
centricities are above ~ 0.02 (estimated from Equation [HJ, the 
extent of the clearing will be greater than this because of the ex- 
tended chaotic zone. Hence, i f the disc partic les have low eccen- 
tricities, the mass estimated bv lSu et alj j2009l) will be correct, and 
the planet mass will be around 10 Jovian masses. However, if the 
disc particles' eccentricities are higher than ~ 0.02, the same clear- 
ing will be achieved with a smaller planet mass. We quantify this 
effect in Figure|5] This compares the shape and extent of the clear- 
ing due to planets of 2, 4, 6, 8 and 10 Jovian masses, under the 
assumptions that particles are removed only in the classical chaotic 
zone (dotted lines) or also in the extended chaotic zone (solid line). 
In each case, particle eccentricities were distributed uniformly be- 
tween and 0.1. We also show as vertical dotted lines the location 
of the edge for particles on circular orbits. We see that the profiles 
for 8 and 10 Jovian mass planets under the Wisdom prescription 
(whose medians straddle 90 AU) are similar to the profiles for the 
4 and 8 Jovian mass planets, respectively, under the extended zone 
prescription. Hence, the planet masses required to achieve a given 
clearing may be less than predicted by the Wisdom criterion by as 
much as 50%. 

Thus, if the edge of the outer disc is located at 90 AU, and the 
particles in the disc have eccentricities of order 0.1, then the mass 
of planet b would be in the range 4-8 Jovian masses, rather than the 
8-10 Jovian masses that it would be if the particles' eccentricities 
were below 0.02. For now it is not possible to make such firm con- 
clusions, since observations have yet to resolve the inner edge of 
the disc, and its 90AU location inferred from SED modelling is de- 
generate with the assumptions made about the particle properties. 
Nevertheless, the inner edge will likely be resolved by future ob- 
servations, and this will help to reduce t he uncertainty in th e planet 
mass. For example, some authors, e.g., ICurrie et al.l ( |20I1|) . claim 
a lower mass for planet b of 6-7 Jovian masses, and this is consis- 
tent with the planet truncating a disc of planetesimals, whose ec- 
centricities range up to 0.1, at 90 AU, if the edge indeed be located 
there. The eccentricities of the disc particles may have been excited 
to such an extent by the g rowth of large planetesimals in the disc 
jKenvon & Bromlevll201Cll) . secular perturbations fro m the planets, 
if one or more is eccentric jMustill & Wvatt|[2009h . or sweeping 
by mea n motion or secular resonances during past evoluti on of the 
system jMustill & Wvatll201 itlMinton & Malhotr j|201 1[) . 

Finally, we note that our main conclusion - that the mass of 
planet b estimated from the extent of the chaotic zone will depend 
on the eccentricities of the disc particles - does not depend on the 
actual location of the disc edge. If the disc prove to be at a greater 
radius, the mass estimated for a disc of eccentric particles will still 
be greater than that estimated for particles on circular orbits, al- 
though each will be greater than its corresponding mass estimated 
for a disc edge at 90 AU. 

4.3 Post-Main Sequence evolution of planetary systems 

There are also important implications for the evolution of planetary 
systems when a star loses mass during post-Main Sequence 
evolution. When stars lose mass on the Asymptotic Giant Branch 
the planet: star mass rat ios increase and prev i ously s table systems 
can be destabilised (Debes & SigurdssonI |2002| : iBonsor et al.l 
I20T ij). Stars typically lose mass at relatively modest rates, and 
under these conditions orbits ex pand adiabatically and the ratios of 
semi-major axes are unchanged dVeras et al.l20T l]). However, since 
the planet:star mass ratio increases as the star loses mass, quantities 
which depend on this ratio, such as the Hill's radius and the widths 
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Figure 5. Comparison of disc edge profiles for the disc around HR 8799 
truncated by planet b at 68 AU. The maximum eccentricity of the disc par- 
ticles is 0.1. Solid lines show disc profiles derived from the encounter map 
results; dotted show profiles when truncated by the classical Wisdom pre- 
scription. The vertical dotted lines show the location of the disc edge for 
particles on perfectly circular orbits. Planet masses increase from left to 
right: 2, 4, 6, 8 and 10 Jovian masses. Note that all the disc surface densities 
are 10* far from the planet; the different heights visible here ai'e due to the 
effects of the 2:1 resonance at 108 AU. 

of resonances and the chaotic zone, will change too. The chaotic 
zone expands, and after mass loss particles which were previously 
on stable orbits may find themselves in the chaotic zone. The pres- 
ence of this reservoir of newly unstable material may explain the 
existence of metal pollution in the atmospheres of White Dwarfs 
(Zuckerman et al. 2003), as well as hot discs orbiting some White 
Dwarfs jFarihi et alj |2009[) . as bodies a re scattered onto highly 
eccentric orbits a n d tidally dis r upted jZuckerman et alj |2003| ; 
iBonsor et alj|201 ll) . iBonsor et al.l ( I2OIII) estimated the amount of 
destabilised material simply by expanding the /^^'^^ zone according 
to the new, lower, stellar mass. While this will correctly describe 
the new chaotic zone at low eccentricity, seemingly vulnerable 
particles with higher eccentricity may escape being engulfed by 
the new chaotic zone. This is because the lower boundary of the 
extended chaotic zone increases as ii increases . Howe ver, a further 
complication may counteract this: IVeras et ^ j201 ih showed that 
if mass loss rates are relatively rapid then particles' orbits do not 
expand adiabatically and they may acquire some eccentricity. This 
could push seemingly safe bodies into the extended chaotic zone. 
Since the eccentricities required to enter the extended chaotic zone 
are relatively modest, this could be an important mechanism for 
destabilising bodies. 
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